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Abstract. The hyperspherical adiabatic expansion is combined with complex scaling and used to calculate 
low-lying nuclear resonances of 12 C in the 3a-model. We use Ali-Bodmer potentials and compare results 
for other potentials a — a with similar 8 Be-properties. A three-body potential is used to adjust the 12 C- 
resonance positions to desired values extending the applicability of the method to many-body systems 
decaying into three a-particles. For natural choices of three-body potentials we find 14 resonances below 
the proton separation threshold, i.e. two + , three 2 + , two 4 + , one of each of 1*, 2", 3*, 4", and 6+. The 
partial wave decomposition of each resonance is calculated as function of hyperradius. Strong variation 
is found from small to large distance. Connection to previous experimental and theoretical results are 
discussed and agreements as well as disagreements are emphasized. 

PACS. 21.45. +v Few-body systems - 21.60.Gx Cluster models - 25.70.Ef Resonances - 27.20. +n. 6 < = 

A <= 19 



1 Introduction 

The low-lying nuclear bound states and resonances of 12 C 
have been the subject of numerous investigations. Surpris- 
ingly a number of issues are still not settled, e.g. what are 
1 the energies of the low-lying resonances below the pro- 
, ton separation threshold 15.96 MeV of excitation energy, 
what are their angular momenta, their structure in gen- 
eral, and their decay properties. Substantial experimental 
efforts are presently devoted to find answers to these ques- 
tions [1] . This is partly motivated by the interest from as- 
trophysics [2] , but the research has basic intrinsic interest 
in its own right. 

Experiments often focus on specific properties and/or 
specific resonances, but taken together more and more ac- 
curate data accumulate [3] . Open questions remain on + 
and 2 + resonances [4] above the first + resonance at 7.6 
MeV, e.g. in [3] there are only tentatively assigned can- 
didates. Recently this has been addressed experimentally 
both by inelastic scattering of alpha-particles on 12 C [5, 
6], and by studies of the beta-decays of 12 N and 12 B [7]. 
These studies are in reasonable agreement on the existence 
of a broad + state near 10 MeV, but give three different 
suggestions for the position and width of the first 2 + res- 
onance. 

The accurate measurements and the related careful 
analysis call for a reliable theoretical description. The the- 
oretical efforts go back at least fifty years but consensus 
has not been reached. Advanced recent attempts can be 
divided into microscopic models dealing either directly 



with an unbiased 12-body nuclear problem or with vari- 
ous cluster constraints of four nucleons within each of the 
three a-particles. Within the first group we find the no- 
core shell model [8] , the stochastic variational method [9] 
and Greens function Monte Carlo calculations [10]. None 
of these methods have so far been able to provide satis- 
factory answers to the present continuum problem. In the 
second group we find antisymmetric molecular dynamics 
[11], Fermionic molecular dynamics [12], and resonating 
group methods [13]. These methods are very ambitious 
as well but more directed to account for possible cluster 
structures. However, they have also not provided convinc- 
ing answers as for example evidenced by mutual disagree- 
ments. Furthermore, the same tendency of concentrating 
on specific properties and resonances is seen as in the ex- 
perimental studies. Systematic investigations of all low- 
lying resonances with a given model are very few [15,16, 
14,11], and not easily reconciled with the newest experi- 
mental information. 

These microscopic theories have so far not provided en- 
ergy distributions of the a-particles after decay of the res- 
onances. These observables are difficult to compute with 
a reasonable accuracy, but nevertheless this is precisely 
where the accurate and complete experimental informa- 
tion accumulates [7]. To address these problems we scale 
down the ambitions to a practical yet challenging size 
where systematic investigations are possible. —3a cluster 
models suggest themselves since energies below the proton 
separation threshold at 15.96 MeV only allow 3a and 7- 
emission. Many-body resonances decaying into three par- 
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tides necessarily reduce to a three-body problem at large 
distance. The validity of 3a models extends for short- 
range interactions to surprisingly small distances where 
the three particles almost touch each other [17-19]. The 
a-cluster structure may dominate and the description can 
possibly be extended to even smaller distances. However, 
detailed information at small distance can not be expected 
without introduction of nucleonic degrees of freedom. For 
this reason three— a models can only be expected to de- 
scribe intermediate and large-distance properties related 
to a-clustering. In particular electromagnetic transitions 
depend sensitively on short-distance properties of the wave- 
functions and are therefore in general not reliably acces- 
sible within these models. 

The limitations of these models are not established. 
The energy distributions after decay are large-distance 
properties and as such possible to compute rather accu- 
rately with few-body techniques [20]. However, even this 
3a problem is very challenging. First because we are then 
dealing with a three-body quantum mechanical Coulomb 
problem in the continuum. Second because we want a 
comprehensive picture describing all resonances and all 
three-body observables within the same model. The sec- 
ond requirement is especially demanding since the res- 
onance may vary from two-body to genuine three-body 
large-distance structures. Coulomb and short-range inter- 
actions must then be treated on an equal footing. However, 
the techniques and methods are available and recently ap- 
plied in similar practical three-body computations [21]. 

We shall use the hyperspherical adiabatic expansion 
method [22] combined with complex scaling of the coordi- 
nates [19]. This method restricts the applications to res- 
onances where the width is relatively small compared to 
the real part of the energy. In any case as these states are 
closest to the real axis they have the strongest direct influ- 
ence on observables. Furthermore, as obtained by analyt- 
ical continuation they have the smallest uncertainty due 
to the choice of parametrization of the interactions. To 
describe properly the different types of asymptotic large- 
distance structures we use the Faddeev decomposition. 
Then the individual components can simultaneously ac- 
count for several two-body substructures like narrow reso- 
nances in corresponding subsystems. This is a tremendous 
advantage over one-component methods. 

The purpose of the present work is to (i) establish that 
the 3a-cluster model can answer questions related to prop- 
erties of many-body resonances, (ii) give a survey of the 
structures of possible 12 C-resonances below the excitation 
energy of 15.96 MeV, and (hi) lay out the foundation for 
calculations of energy distributions after decay. We shall 
build on experience gained from three-body investigations 
of bound states and continuum properties of 6 He [23], n Li 
[25,24,20], 17 Nc [26], 6 Be and 6 Li [27] as well as more 
general investigation of three-body resonance properties 
[28-30]. In section 2 we shall first give a brief sketch of 
the theoretical formulation. In section 3 we discuss the 
resonances and their structures for a specific two-body 
interaction. In section 4 we describe the sensitivity to in- 



teractions. Finally section 5 contains a summary and the 
conclusions. 



2 Basic theoretical ingredients 

We use techniques described in details in previous publica- 
tions. It suffices here to give a brief sketch of the procedure 
employed in the present work. The general framework can 
be found in [22]. The new results found for 12 C will be 
discussed in more details in the following sections. The 
adiabatic hyperspherical expansion method with the Fad- 
deev decomposition is well established in applications to 
nuclear three-body systems [25]. The combination with 
complex rotation to compute resonances is also known to 
be very efficient for these systems [19]. The method con- 
sists of a number of steps: First we define the coordinates 
where the hyperradius p is the most important, i.e. 

p 2 = \jZ^-^f-^jl^-Rf > (!) 

i<j i—1 

where r j is the coordinate of the a-particle number i and 
R is the center-of-mass coordinate of 12 C. The factor of 4 
is arbitrary and chosen to correspond to a normalization 
mass equal to that of the nucleon. The remaining relative 
coordinates are all dimensionless angles. 

Second we choose a two-body interaction as input for 
the angular part of the complex rotated Faddeev equa- 
tions. This interaction should reproduce the low-energy 
scattering properties of all pairs of particles in the three- 
body system. A smaller amount of data like scattering 
length and effective range may also be sufficient, or per- 
haps the low-lying two-body resonance energies and their 
widths. We use the Ali-Bodmer potentials [31] in various 
combinations, and sometimes compared with results from 
even simpler interactions adjusted to reproduce selected 
8 Be-propcrties. The hyperradius is fixed and the solutions 
to the angular wavefunction <P n are calculated. 

Third a three-body Gaussian potential, V36 = S*exp(— p 2 
is selected. Such potentials are necessary to reproduce ac- 
curately few-nucleon observables [32] . We choose b « 6 f m 
such that p = b corresponds to three touching a-particles. 
This three-body potential is diagonal, i.e. added to each of 
the adiabatic potentials whereas we assume that the cor- 
responding three-body couplings are zero. This construc- 
tion maintains the structure of a three-body state, but by 
varying the strength the energy position can be adjusted 
to reproduce the measured value for each resonance. 

An individual adjustment is crucial in a comparison of 
partial widths which depend exponentially on the energy 
due to barrier penetration. Since the small-distance struc- 
ture easily can vary from state to state the three-body 
potential should in principle be state-dependent. However, 
we choose the same potential for all states of given angu- 
lar momentum and parity, J 71 '. One simple reason is that 
this essentially corresponds to the number of known ex- 
perimental data and we can then estimate or predict the 
positions and widths of other states with the same J w . 
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States with different J 71 " can be expected to differ more 
in structure and thus require different three-body poten- 
tials. We choose to ignore the two bound states because 
they, and especially the ground state, have a smaller size 
and their structure can then be expected to involve more 
nucleonic (less a-cluster) degrees of freedom. Although it 
is possible to construct a three-body potential which re- 
produces energies of both the ground state and the lowest 
+ resonance [33,34] it is uninteresting in the present con- 
text since the important intermediate and large distance- 
structure would remain unchanged. 

Fourth we find the radial wavefunctions, f n {p), by solv- 
ing the coupled set of radial equations arising from expan- 
sion of the total wavefunction & on (p n , i.e. 

^^iE^'W)' ( 2 ) 

" n 

The expansion coefficients f n (p) are exponentially decay- 
ing for resonances when the rotation angle 8 of the hy- 
perradius is larger than corresponding to the three-body 
resonance. Then both real and imaginary parts, E = 
Er — iEj 7 of the resonance energy, Eq, are determined 

by f n (p — * oo ) = C n exp(+zKp) with k = \J IttiEq / h 2 . 

In principle some of the channels could correspond 
to two-body bound states or narrow resonances which 
asymptotically would have radial wavefunctions decaying 
with a wave number corresponding to the three-body en- 
ergy minus that tied up in the two-body system. In prac- 
tice this only has marginal importance since the radial 
wavefunction in any case decreases towards zero for large 
distances. The interpretation, and maybe the analysis, in 
terms of direct or sequential decay may however be very 
different. 

Fifth the structures of the resonance wavefunctions 
are computed and expressed in terms of two-body par- 
tial waves. We extract the amplitudes (above called the 
radial wavefunctions) for each of the adiabatic potentials, 
and then we partial wave decompose each of these adia- 
batic components. The dependence on hyperradius can be 
substantial corresponding to a dynamical evolution of the 
resonance from small to large distances. 

In these steps the problems with the Coulomb interac- 
tion are not mentioned although the large-distance asymp- 
totic behavior is mathematically unknown for continuum 
states. We adopt the pragmatic procedure to treat the 
Coulomb interaction completely numerically. The general 
behavior and convergence properties can then be inves- 
tigated and extracted if necessary [21]. Problems would 
reveal themselves intrinsically in the process as numerical 
inconsistencies. 



g-waves. The 8 Be resonances of + , 2 + , 4 + have ener- 
gies (E R ) and widths (r = -2£V) given by {E R ,r) = 
(100,0.010) keV, (2.7,1.5) MeV, and (9.7,8.0) MeV, re- 
spectively. The definitions correspond to poles of the S- 
matrix, and should be compared to the experimental val- 
ues of (E R ,r) = (92,0.0068) keV, (3.1,1.5) MeV, and 
(11.5,4) MeV, respectively. The 4 + -state is not very well 
reproduced although the information is derived from phase 
shifts reproducing the data. 



3.1 Partial waves 

For each of the 14 resonances angular momentum and 
parity constrain the contributing two-body components 
within the three-body system. The partial wave orbital 
angular momenta between two particles in one Jacobi sys- 
tem and that of the third particle are denoted by £ x and £ y , 
respectively. The number of partial waves is rather small 
because the a-particle has zero spin and the symmetry re- 
quirement eliminates in addition several components. The 
Faddeev components are the same in all three Jacobi sys- 
tems. The number of basis states needed to describe the 
necessary number of partial waves can only be decided a 
posteriori. Accuracy is optimized by choosing large basis 
sets, enumerated by the hyperspherical quantum number 
K maxi when the contribution is large for a given partial 
wave. 

Convergence has been achieved with the sets of quan- 
tum numbers shown in table 1. The large values of K max 
should roughly be divided by two to give the correspond- 
ing number of basis states in each partial wave. This means 
that the total number of basis states is about 3 x 90 for 
£x = £y = 0, and in total for the + -resonanees the basis 
consists of about 700 states. Some of the other resonances 
have a larger number of basis states. With present day 
computers these numbers are far from being large but still 
extremely efficient compared to the use of only one Jacobi 
system and all possible partial waves consistent with one 
value of K max [35]. These high partial waves are needed 
for one Jacobi system to describe configurations with two 
spatially close-lying particles relatively far from the third 
particle. The Faddeev decomposition takes care of that in 
the present formulation. 

A given angular momentum and parity of an intended 
12 C-state needs a set of these partial wave components. 
However, more than one 12 C-state with this spin and par- 
ity may be found from such a set. The probability to find 
each partial wave for each resonance is also given in table 
1. We shall postpone the discussion of these values until 
section 3.4. 



3 Resonances and their structure 

The sequence of resonances is computed with the Ali- 
Bodmer potential "AB(a')" [31] which reproduce the s, 
d and qualitatively also the g-w&ve a-a phase shifts. The 
interactions for higher partial waves are the same as for 



3.2 Adiabatic potentials 

The adiabatic potentials result from a full quantum me- 
chanical solutions of the angular part of the Faddeev equa- 
tion for fixed hyperradius. Each potential corresponds to 
a specific combination of the partial waves, in some cases 
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Table 1. Components included for each J w state of C. The 
columns two and three give orbital angular momenta in each 
Jacobi system column four gives the maximum value of the hy- 
permomentum K. In the last columns Wi give the probabilities 
in % for finding these component in the i-th resonance. 
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strongly varying with p. These structures form the ba- 
sis for the radial solution which truly is a resonance in 
the three-body continuum. Very few angular eigenvalues 
are usually needed for convergence, and their behavior are 
decisive for the resulting resonances. This is not a trivial 
conclusion, since each point of these potentials by defini- 
tion corresponds to the same p, but different configura- 
tions are otherwise allowed. It is not obvious at all that 
the structures of the resonances predominantly arise from 
these combinations of configurations contributing to each 
adiabatic potential. 

We show in fig. 1 the real parts of these potentials in- 
cluding Vj,b corresponding to the quantum numbers of each 
of the resonances. The imaginary parts arc small, oscillate 
around zero and vanish at large p. All these potentials di- 



verge as p~ 2 when p — > due to the generalized centrifugal 
barrier. They vanish as 1/p for large p where Coulomb is 
the only contributing interaction. At relatively small dis- 
tances the potentials have minima supporting the bound 
states and resonances. The sizes of these pockets are larger 
than obtained from the two-body interaction which in 
most cases is far from sufficient to place the resonance 
at the correct energy. Therefore we added in each case 
a diagonal, short-range three-body potential designed to 
mock up effects of the intrinsic degrees of freedom in a 3a- 
model at short distances where the three a-particles over- 
lap. This allows us to stay within the three-body model 
and still provide the proper boundary condition at small 
distance at the right energy. 

The three-body potentials are individually adjusted to 
reproduce one of the energies, but chosen to be the same 
for all adiabatic potentials of given angular momentum 
and parity. We typically choose the best known of the res- 
onances since it is the continuum state of lowest energy, i.e. 
closest to the threshold for 3 free a-particles. Obviously 
the bound states are better known but their structures 
can easily differ even more from the 3a-structure than the 
resonance states. The three-body potentials could then 
differ substantially and predictions for other resonances 
probably would be less accurate. The systems with wave- 
functions localized outside the overlap region for the a- 
particles should need less three-body potential and then 
be more reliably calculated in the 3a-model. For exam- 
ple the ground state could be a spatially confined struc- 
ture while + -resonances could be fairly good three-body 
states. 

For the + , 2 + and 4 + -cases two or three minima ap- 
pear in different potentials indicating the possibility of dif- 
ferent structures of low-lying resonances. For the l + -case 
the minimum only appears after addition of the three- 
body potential indicating that the corresponding struc- 
ture at small distance differs substantially from three a- 
particles. All negative parity states exhibit one minimum 
at distances larger than the 6 fm chosen as the range of 
the three-body interaction. 

3.3 Resonance energies 

The length coordinate, p, is rotated in the complex plane 
and the adiabatic potentials including non-diagonal terms 
are computed. By adding the diagonal three-body poten- 
tial all quantities are specified in the coupled set of radial 
equations. The solutions are the three-body resonances 
where key quantities are real and imaginary parts of the 
energies. The strength of the three-body potential is used 
to move one bound state or one resonance into a preferred 
position for each set of angular momentum and parity. If 
there is more than one solution they are moved simultane- 
ously by the same potential. The range of the three-body 
potential is kept at the same value implying that the width 
is predicted and not fitted. The associated uncertainties 
are discussed in section 4. 

The computed results for the resonance energies are 
shown in fig. 2. The individual figures all exhibit the same 
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Table 2. Calculated and measured energies En (in MeV) 
and partial a-decay widths -Tr (in keV) of the resonances for 
different J*. Since our formalism cannot describe an isospin 
one state, we found only one 1 + resonance whose energy has 
been adjusted to reproduce both measured values. Experimen- 
tal values (labeled "exp") are from [3,7,1] and calculated re- 
sults (labeled "th") are obtained with the two-body interac- 
tion "AB(a')" and the three-body interaction parameters S 
(in MeV) and b (in fm) given here. The energies are measured 
from the 3a threshold. The computations are for a rotation 
angle 9 (in rad). The 1+ state at 7.84 MeV labeled T = 1 has 
isospin one. 
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dominant features with a series of points forming almost a 
straight line. Each of these points corresponds to a contin- 
uum state. These discrctized three-body continuum states 
appear due to an imposed condition of zero radial wave- 
function at the boundary of a large box. The energies ap- 
pear on the straight line arising by rotation of an angle 29 
around the origin. The fluctuations are due to numerical 
uncertainties. 

The points above the straight line are the resonance en- 
ergies appearing with positive real parts, Er, and negative 
imaginary parts, Ei = —T/2, where r is the width of the 
resonance. The radial wavefunctions of these resonances 
vanish exponentially at large distances corresponding to 
outgoing waves in all channels. This boundary condition is 
identical to those valid for the bound states which appear 
on the negative side of the real axis. 

The computed resonance energies are collected in table 
2 where the experimental values also are given. Three Up- 
states appear, i.e. the ground state and two resonances. 
The first resonance is very close to the threshold with an 
extremely small width of about 62 eV. It is then difficult to 
distinguish from a discretized continuum state without in- 
formation from its wavefunction. The second + -resonance 
appears around 4 MeV with a width of about 1 MeV. The 
real energy disagrees with previous experiments [3] , but it 
is consistent with the analysis of the recent experiment [7, 
1] as quoted in table 2. 



We find four 2 + -states, i.e. one excited bound state and 
three resonances. Only two resonances are suggested by 
experiments [3,7,1]. With the natural choice of b = 6 fm 
the highest of the computed resonances can easily be ad- 
justed to the experimental value [7,1]. It would perhaps 
be natural to adjust the lowest resonance to the measured 
energy. However, an attempt to do this failed because a 
smaller attraction than given in table 2 move the energy 
upwards but at the same time the top of the confining bar- 
rier is approached and exceeded leading to much faster 
increase of the width. The result is that the width in- 
creases quickly whereas the energy remains essentially un- 
changed. The same behavior is seen for the second of the 
2 + -resonances. Only the third follows naturally the exper- 
imental specifications with this choice of one three-body 
interaction of the present Gaussian form. 

We find two 4+-resonances. If the three-body strength 
is reduced from the value in table 2 the widths of both res- 
onances increase dramatically around 4 MeV and 9 MeV, 
respectively. Therefore the most consistent solution is that 
the highest 4 + -resonance is found experimentally while 
the lowest has been hidden behind other states in previ- 
ous experiments. Consequences of allowing variation of b 
are discussed in section 4. 

Two 1 + resonances are known experimentally, where 
the 7.84 MeV state is an isospin T — 1 state. The two cal- 
culated l + -resonances are artificial in the sense that only 
one of these resonances is found in the calculation, since 
a state with T = 1 cannot be described by our formalism. 
The energy is then adjusted through the three-body po- 
tential to either of the measured values corresponding to 
isospin zero or one. The fact that only one state is found 
in the computation reflects that we only have the correct 
degrees of freedom to describe T — states. The reason 
for also adjusting to the T = 1 energy is the expectation 
that the decay into three a-particles would occur through 
the same intermediate and large-distance structures as for 
T = 0. We do not imply that the T = 1 state has the same 
structure at small distances. 

All other states than = + ,2+,4+ only appear 
once. Their energies are then adjusted by use of the three- 
body potential and the related width for the range of 
b = 6 fm can be compared to the measured values. Since 
only one 2~ and one 4~-resonance appear we suggest that 
the uncertain assignment of the second experimental 2~- 
resonance should be changed to 4~ . The 3 + and 6 + states 
are not known experimentally and we can only make rough 
guesses of the three-body strength. Since 6 + belongs to the 
family of even angular momenta and even parity we used a 
similar strength of —20 MeV. The unnatural parity of 3 + 
puts it in the same class as the 1 + state but the strength 
of —20 MeV already leads to an energy of about 6 MeV. 
This 3 + estimate is then very uncertain. 

It is remarkable that the ratio between theoretical and 
experimental values for the widths in almost all cases is 
around 2. The widths are in most cases fairly well esti- 
mated by the tunneling probabilities through the respec- 
tive barriers shown in fig. 1, see [36]. This corresponds 
to an exponential dependence on energy and deviations 
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within a factor of two can then be considered to be very 
accurate [36]. These results are obtained without special 
treatment of the small distance properties beyond that 
corresponding to the 3a-structures. The possibility of pre- 
formation factors or a-particle spectroscopic factors, as 
appearing in ordinary a-decay, seems to play a minor role. 
It is interesting to note that we systematically calculate 
slightly too large widths indicating preformation factors 
at the surface of about 1/2. 

Two exceptions from the general rule are seen in the 
second + and the two 1+ resonances, where the com- 
puted width is 3.4 times smaller for + and larger for the 
two 1+ states by 2.6 • 10 3 and 5.3 • 10 5 , respectively. The 
1 + states require in contrast to all the other states pre- 
formation (spectroscopic) factors three and five orders of 
magnitude smaller than for the other states. This reflects 
totally different structures corresponding to a full break 
down of the 3ai-descriptions at small distances. This is 
particularly pronounced for the T = 1 state which cannot 
be described by the T = a-particle building blocks. The 
preformation factor correspondingly comes out as exceed- 
ingly small. 

The different structure is less obvious for the lowest 
1 + state with T = which is not forbidden for symme- 
try reasons in the 3a-models. However, a strong hint is 
found in the model by the unusually large lowest possi- 
ble hypermomentum of 8. The immediate guess of K = 1 
for 1 + states are ruled out, together with all other hyper- 
momcnts below K = 8, by the imposed boson symmetry 
of the three a-particles [37]. This is in clear contrast to 
all the other J 71 ' states where the lowest hypermomentum 
consistent with angular momentum always contribute, see 
table 1. 

These results are also, apart form the 1 + anomalies, re- 
markable for at least one other reason, i.e. a monotonic de- 
pendence on energy, angular momentum and parity docs 
not exist. The barriers in hyperradius provided by the 
calculation must then be physically reasonable, in turn 
implying that the concept of the hyperspheric adiabatic 
expansion is in rather close agreement with the physical 
process. The systematic correlation with the experimen- 
tal widths, varying by orders of magnitude, is otherwise 
totally unexplained. It should here be emphasized that 
many of the resonances are not well described as cluster 
states [11]. This is particularly clear for the two 1 + states. 

The three-body potentials with b = 6 fm seem to fall in 
two groups, i.e. one with a strength S of about —20 MeV 
(0+, 2+, 3+, 4+, 6+) and one with S w -2 MeV (1", 2", 
3~, 4~). This should indicate that the amount of a-cluster 
states within each group is roughly the same as also re- 
flected by the comparison of calculated and known widths. 
Exceptions are the 1 + states with much larger strengths of 
about —80 MeV. The computed widths for these states in- 
crease substantially with the range b of the three-body pa- 
rameter when the strengths correspondingly are increased 
to maintain the energy. The apparent spectroscopic factor 
for the isospin zero state would for b « 4 fm be comparable 
to the values for the other states as reported in [36]. 



Beside the energy the size is a characteristic quantity 
for any system. This is straightforward for bound states 
but resonances decay and are therefore infinitely large. 
However, the complex scaling provides a simple measure 
of the extension of the small distance part of the wave- 
function, i.e. expectation values related to the hyperra- 
dius. The usual relation is [22] 

rims = ^\P 2 W) cxp(^) + R 2 a , (3) 

where R a — 1.47 fm and r rms are the root mean square 
radii of the a-particle and 12 C in the 3a-model, respec- 
tively. The finite size of the three a-particles are then ac- 
counted for. Complex rotation by 6 means that the square 
of the hyperradius is multiplied by exp(i29). The expec- 
tation value in eq.(3) is then expressed in the rotated sys- 
tem, and now this quantity is finite. They are still complex 
numbers and the absolute square can be used as the mea- 
sure of size. Another tempting possibility is to interpret 
the real part as the radius of the resonance and the imag- 
inary part as the variation (uncertainty) of the radius of 
the continuum wavefunctions as energy changes through 
the interval Er ± T/2 around the central energy [38]. 

The numerical results for the resonances and bound 
states are given in table 3 for both prescriptions. The real 
parts are in all cases very close to each other and certainly 
within the "uncertainty" measured by the imaginary part. 
For the two bound states the imaginary part reflects the 
small uncertainty in the numerical computation, since here 
the values should strictly be real. The ground state re- 
sult is within the uncertainty identical to the root mean 
square radius of 12 C, i.e. 2.468 fm. The "uncertainties" 
are small for narrow resonances and relatively large for 
broad resonances. The bound 2 + -state is smaller than the 
ground state and all the resonances are larger. The + 
and 2 + -resonances are especially large because they are 
built of the same partial waves and orthogonality then 
apparently push them outwards. The comparison to the 
old results derived from the resonating group method [39] 
and the ones from the generator coordinate calculation 
[16], as well as the recent computation [40] and the anti- 
symmetrized molecular dynamics results from [11] are in 
agreement within the uncertainties. It should be empha- 
sized that a detailed comparison should take the energy 
differences into account. From halo physics of neutral par- 
ticles we know that energy and radius are closely corre- 
lated. The radial variation with binding energy is much 
smaller for charged particles and finite angular momen- 
tum. 



3.4 Resonance structures 

The resonance structure is determined by combining the 
amplitude obtained as the solution to the radial wave- 
function with the individual structure of each adiabatic 
component. The real parts of the radial wavefunctions arc 
shown in figs. 3 and 4 for the two bound states and each 
of the 14 resonances. At least two general features are 
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Table 3. Sizes expressed as root mean square radii (in fm) of the 12 C bound states and resonances. The second column is 
prms = | < p 2 > I 1 / 2 , and the fourth column is without absolute squares of the expectation values. The third and fifth columns 
are root mean square radii of 12 C states obtained by eq. (3) with the two different < p 2 > prescriptions, respectively. The last 
four columns are computations from [11,39,16,40]. The root mean square charge radius of 12 C ground state is < r 2 y 1 / 2 — 
2.4829 ± 0.0019 fm [41]. 
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striking. First, only very few of these wavefunctions con- 
tribute significantly Furthermore, for all states one adi- 
abatic component is much larger than all other contri- 
butions. This reflects the rather fast convergence of the 
adiabatic expansion where three eigenvalues often provide 
an accurate solution. Second, the wavefunctions oscillate 
with exponentially decreasing amplitudes as functions of 
p. This behavior is qualitatively the same for bound states 
and resonances which precisely is the reason for using com- 
plex scaling to compute resonance wavefunctions and en- 
ergies. The imaginary parts are not shown to avoid clutter- 
ing the figures. They exhibit similar oscillatory behavior. 

The probabilities for finding the different adiabatic 
components in each resonance wavefunction are given in 
table 4. The structures of the resonances are then essen- 
tially contained in these dominating angular wavefunc- 
tions, where each is related to one adiabatic potential. For 
each resonance we therefore take the dominating contri- 
bution and decompose into the partial waves specified in 
table 1. These partial wave decompositions are shown in 
figs. 5 and 6. The overall feature is that the dependence on 
hyperradius is rather strong. However, we should empha- 
size that the small distance structure only has little physi- 
cal significance. Numerically the variation arises from cou- 
plings of low- lying adiabatic potentials which change rapidly 
in this region due to the different behavior of s, d and g 
two-body potentials imposed by the different behavior of 
the phase shifts. 

The two dominating wavefunctions for the + -resonances 
have very different behavior. The first is dominated by 
£ x = £ y = at about 10 fm, but at larger p the higher par- 
tial waves all become almost equally important. The large- 
distance behavior is consistent with the description of one 
a-particle far away from a spatially confined 8 Bc(0 + )- 



Table 4. Contribution W„ of each adiabatic potential, labeled 
by n, to the total norm of the wave function of each 12 C-state. 
Only those contributions of more than 1% are given. 
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structure. The higher partial waves arise from s-waves in 
the "unnatural" Jacobi systems. The second wavefunc- 
tion has about equal contributions at small p from £ x = 
£ y = 0,2, but at large p the s-wave component increases 
strongly while the higher partial waves correspondingly 
decrease towards zero. This large-distance structure is con- 
sistent with all three a-particles symmetrically distributed 
far from each other. 

The first of the other natural parity wavefunctions, 
corresponding to 1~, 3~, 2 + , 4 + and 6 + , all have £ x = 
as the dominating component at large distance. This is 
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consistent with a 8 Bc(0 + )-structure and the third particle 
with l y equal to the total angular momentum. Several 
other components are small and of comparable magnitude. 
The lowest + -wavefunction also falls into this category. 

The first wavcfunctions for unnatural parity, 1 + , 3 + , 
2~, AT , all have £ x = 2 as the dominating component both 
at small and large distances. The angular momentum, £ y , 
of the third particle for this component is as small as pos- 
sible consistent with the total spin and parity. Only one 
more component gives significant contributions for each of 
these states, i.e. £ x = 4 combined with the lowest possible 
l y for 1+ and 4", and t x = 2 and £ y = 3 for 2". The latter 
component is larger than the contribution from the higher 
partial wave of £ x — 4 and the lowest possible Rvalue of 
3. 

The third 2 + -wavefunction is similar in structure with 
one dominating component of £ x = 2 and £ y = 0, and an- 
other component of half the size with £ x = 2 and £ y = 4. 
Here the component with £ x — 2 and £ y = 2 only con- 
tributes insignificantly. The second 2 + -wavcfunction has 
three about equal contributions at large distance from 
(£ x ,£ y ) = (0,2), (2,0), (2,2). The second 4+-wavefunction 
also has three comparable contributions from (£ x ,£ y ) = 
(0, 4), (4, 0), (2, 2), and in addition several smaller compo- 
nents. 

We emphasize that the contributions from the indi- 
vidual adiabatic potentials should be combined with their 
respective amplitudes to give the total structure of the 
resonances. These amplitudes depend on hyperradius and 
the probabilities in table 4 are averages over all p. The ex- 
ponential decrease with p then obviously heavily enhances 
the structures at the small distances. The relative distri- 
butions at larger distances may change substantially. 

The different partial waves receive contributions from 
different adiabatic components. The probabilities given in 
table 2 can be obtained by combining the information in 
figs. 5 and 6, and figs. 3 and 4. The results reflect the 
conclusions from the above discussion. 



4 Interaction dependence 

We have to distinguish between two types of inaccuracies, 
i.e. arising from (i) inaccuracies in the numerical computa- 
tions for given sets of input parameters, and (ii) model as- 
sumptions and uncertainties in the corresponding param- 
eters. Here the uncertainties in (i) originate from cut-offs 
in different basis expansions, i.e. the partial waves and the 
number of Jacobi polynomials for each of the three Fad- 
deev components, and the number of adiabatic potentials. 
The convergence can be seen directly in the decreasing am- 
plitudes for both partial waves and adiabatic potentials. 
We emphasize that the contributing terms vary with hy- 
perradius as discussed in the previous section. Also the ro- 
tation angle should be chosen as small as possible to opti- 
mize accuracy while still allowing distinction of the three- 
body resonance from the background continuum states. 

Model uncertainties are more difficult to assess. We 
shall first consider those intrinsic to the cluster model, 



and afterwords compare to experiments and other more 
microscopic computations. 

4.1 The 3a-cluster model 

The model uncertainties are essentially related to the dif- 
ferent choices of interactions reproducing the low-energy 
a — a scattering properties. The Ali-Bodmer potentials 
with the correct number of two-body bound states are 
constructed with repulsive cores to simulate the effect of 
Pauli forbidden states. The radial dependence then varies 
from being strongly repulsive to fairly attractive implying 
rapidly varying adiabatic potentials with (avoided) cross- 
ings. To test the sensitivity we constructed a purely at- 
tractive potential "attr" with the two-body complex reso- 
nances energies, (E R , r) = (100, 0.01) keV, (3.2, 1.4) MeV, 
and (11.6,3.9) MeV in £ = 0,2,4 waves, respectively. 
These values are rather close to the measured values. We 
furthermore constructed the potential "AB(a'd')" by com- 
bining the two-body s and d-waves from another Ali-Bodmer 
potential while substituting the g-wave by the above at- 
tractive potential. This potential has £ = 0, 2, 4 8 Be-resonances 
at (E R , r) = (90, 0.006) keV, (3.0, 1.4) MeV, and (11.6, 3.9) 
MeV, respectively. 

The angular eigenvalues for different potentials are 
compared in fig. 7 for some of the resonances. The striking 
feature is that the large-distance behavior of the spectra is 
the same. First the eigenvalues increase linearly with hy- 
perradius simply because the leading order is the Coulomb 
potential multiplied by p 2 . Second the coinciding spec- 
tra at large distances reflect that the resonance positions 
are crucial. These in turn essentially only depend on the 
two-body scattering lengths. Both spectra and resonance 
energies are to leading order then roughly model indepen- 
dent as determined by the low-energy two-body scattering 
properties. 

The complex resonance energies are compared in fig. 
8 for different interactions. First for the + -resonances 
we modestly only change (from AB(a') to AB(a'd')) the 
s-waves a little while both d and ^-potentials are quite 
different. The result is that the energies essentially are 
unchanged for these potentials. The recently found ex- 
perimental result has an unusually large width of about 
3.4 MeV [7,1]. This is somewhat larger but still compara- 
ble with the widths found from inelastic a — a scattering 
experiments [5,6]. We then investigate the dependence on 
the three-body interaction. We follow the resonance as 
function of the three-body strength for fixed values of the 
range. As the energy increases also the width increases 
almost along the continuum background spectrum. If we 
instead keep the energy of the lowest + -resonance by cor- 
related changes of three-body strength and range, the sec- 
ond + -resonance moves relatively little in position while 
the width increases strongly. If both the measured posi- 
tion and the unusually large width should be reproduced 
the range has to be very large and the necessary rotation 
angle would not be numerically accessible. 

For the 2~-case we compare results from the same two 
interactions where we find one resonance at roughly the 
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same position and width in both case. This is more strik- 
ing now where the dominating components have £ x = 2 in 
contrast to + where the similar £ x = O-components domi- 
nate. The variation with three-body strength is rather sim- 
ilar. The resonance position can easily be adjusted with 
reasonable parameters in a fairly large interval around 
5 MeV. Thus even quite different potentials with the same 
two-body resonances give roughly the same three-body 
structure. 

For the 2 + -resonances we compare in fig. 8 the results 
for AB(a') with the results from, attr, a very different po- 
tential with only attraction in all partial waves. First we 
note that the resonance positions vary differently with the 
three-body interaction when we use the initial Ali-Bodmer 
potential. The lowest only moves a little and essentially 
parallel to the background continuum. The second reso- 
nance has a strong tendency to move into the continuum 
reflecting the quickly increasing width as a consequence 
of a decreasing barrier. The third resonance moves faster 
but in the region of the measured energy and width. 

Second it is remarkable that the very different purely 
attractive potential also leads to one bound state and 
three resonances in the same energy region. Again the low- 
est resonance is very close to the 3a-threshold, and even 
below for a fairly weak attractive three-body potential. In 
this case it would seem natural to use a repulsive three- 
body interaction but then the bound state is too easily 
pushed up into the continuum. The level spacing is too 
small as for most other cluster models. The variation with 
three-body potential is similar to the behavior seen for the 
other potential. Now the second and third resonance po- 
sitions are somewhat closer. From these energy variations 
it is then again natural to associate the highest-lying to 
the experimentally found resonance [7,1], and fix its po- 
sition by adjusting the three-body strength for different 
values of the range. The width of this resonance and both 
positions and widths of the other two resonances are then 
determined. 

For the 4 + -resonances and the initial Ali-Bodmer po- 
tential we notice the same behavior as for the 2 + -resonances, 
i.e. the lowest tends to increase the width quickly and dis- 
appear into the background whereas the highest moves 
faster with the three-body strength in the region of the 
measured energy and width. This is again due to the struc- 
ture of the potential. The resonance energy in the mini- 
mum is pushed up towards the top of the barrier where the 
width increases dramatically. The highest 4 + -resonance 
can then, in contrast to the lowest, rather easily be moved 
to the experimental position by adjusting the three-body 
strength. 

For the purely attractive potential we have to choose 
a repulsive three-body interaction if no three-body bound 
states should appear. With the same range of b = 6 fm 
the dominating potential has a high and thick barrier lead- 
ing to a very small width. Instead of an arbitrary adjust- 
ment of the range for this particular case we use the same 
ranges and a similar attractive three-body potential. The 
4+ bound state in fig. 8 should then be considered spu- 
rious and discarded. Then again two resonances appear 



with energies below 8 MeV and widths varying from al- 
most zero to almost 2 MeV. This is within a factor of two 
from the results of the initial Ali-Bodmer potential. 

The structures of the resonances are reflected in their 
partial wave decompositions. The large-distance proper- 
ties of the eigenvalues are almost identical for interactions 
with the same two-body resonances, and consequently the 
partial wave components also coincide for large hyperradii. 
However, it is also obvious that when the short-distance 
properties are qualitatively different then the adiabatic 
potentials must reflect those differences. 

We compare the quantitative differences in fig. 9 for 
some of the adiabatic potentials. For + -statcs the s and d- 
waves dominate and the change into an attractive (/-wave 
has only a small effect except at very small hyperradii. For 
the 2 + -states we compare results where all partial wave 
two-body interactions are different. Nevertheless, the par- 
tial wave decompositions for the three lowest adiabatic 
potentials are remarkably unaffected. The same behavior 
of invariant large-distance properties and relatively little 
variation at small distance is found for the 4 + and 2~- 
states. 

The other resonances are all easily moved to the de- 
sired position by the three-body interaction, and the widths 
are strongly correlated with the measured values, as in 
table 2. The structures of these resonances are also very 
insensitive to the specifics of the interactions. The 6 + - 
resonance seems consistently to fall at an energy of about 
6 MeV which is in the upper end of the considered energy 
interval. This estimate has about the same uncertainty 
as the other natural parity even angular momenta states. 
They are all built on the same low angular momentum par- 
tial waves, see table 1. The additional uncertainty is due 
to the less known two-body interactions of higher angu- 
lar momentum. The 3 + -resonance is much more uncertain 
since the necessary three-body potential could vary from 
almost vanishing to values similar to that of 1 + . 

We can now summarize how resonance energies and 
structures are influenced by two-body interactions repro- 
ducing the low-energy scattering properties. First the num- 
ber of resonances remain unchanged within approximately 
the same region of complex energies. Their energies may 
vary by modest amounts. One exception is the 4+-states 
with the purely attractive potential where the bound state 
has to be discarded as Pauli forbidden. The structures of 
the adiabatic potentials are rather insensitive to the spe- 
cific interaction except at small hyperradii corresponding 
to distances where the three a-particles overlap substan- 
tially. The relative contributions of the adiabatic poten- 
tials may vary at small distances where the contributions 
to the resonances are rather small. 

In particular, the insensitive properties are that (i) the 
bound and the lowest + -statcs are related to the same 
adiabatic potential differing by one radial oscillation, (ii) 
the bound and the second 2 + -states also are related to 
the same adiabatic potential differing by one radial oscil- 
lation, (iii) the second + -resonance and the first (bound) 
2+ state have no radial nodes and both are related to 
adiabatic potentials different from those associated to the 
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lowest resonances, (iv) the third 2 + -state has no radial 
nodes and it is related to a third adiabatic potential, (v) 
all states only appearing once are related to the lowest 
adiabatic potential without radial nodes. 



4.2 Other models 

All angular momenta and parities are also considered in 
some previous model computations. The basic problem in 
comparison is that the interactions often arc different in 
the many investigations. However, the parameters are al- 
ways chosen to reproduce some measured quantities. The 
shortcoming of such an approach is in general that misin- 
terpreted data inadvertently can be used both to fit the 
interactions and afterwords to compare to model results. 
To illustrate the variation of the model predictions we 
discuss some results from selected representative methods 
which all to some degree include the underlying micro- 
scopic nuclear structure. 

The formulation from antisymmetrized molecular dy- 
namics [11] seems first to aim at reproducing the bound 
state properties. The low-energy excited spectrum then 
differs from ours only by one less 2 + -resonance. The level 
spacing is a little larger, and the 1 + and the last 4 + -states 
are found close to 17.5 MeV. We find all our states below 
16 MeV and the spacing between the two bound states too 
small. Our philosophy is that the ground state has contri- 
butions beyond the 3a-cluster structure and we should not 
necessarily reproduce its energy. In the newest preprint 
[11] the spectrum has changed reflecting the use of dif- 
ferent interactions. Now the number of resonances with 
given angular momentum and parity below 16 MeV are 
the same as in the present model. 

The method of fcrmionic molecular dynamics is less 
developed in practical applications [12]. They find one less 
2+, 1+ and 4~-resonance than us below 15 MeV. 

An older publication using the generator coordinate 
method found many resonances below 15 MeV [15,16]. 
They found one less 2 + than us, and otherwise the same 
number of states in this spectrum. The generator coordi- 
nate method has recently been used extensively to inves- 
tigate the structure of light nuclei [13]. The focus is here 
on the bound states and the lowest resonances. The two 
bound states are too close-lying, and the number of reso- 
nances are too few both compared to experiments and to 
our results. Below 15 MeV they find only one + , no 1 + , 
one 2 + , and no 2". Furthermore, their sequence has very 
little resemblance with the measured spectrum. 

The algebraic cluster model is based on group theory 
[14]. Two more states have appeared in the last publi- 
cation close to an excitation energy of 15 MeV. Below 
15 MeV they find the same spectrum as us except for 
one less 4+-state. The overall spacing is therefore almost 
the same. The spectrum is computed from a model with 
3 a-particlcs in an equilateral triangle where rotation- 
vibration interactions, Coriolis forces and vibration- vibrati 
interactions are needed. 



4.3 Experiments 

In confronting the results in table 2 with experimental 
input the predicted low-lying 2 + -statcs and the lowest 
4 + -statc stand out as problematic. However, experimen- 
tally there is not presently a consensus on the position of 
2 + -resonances in 12 C; as previously mentioned different 
probes find states at different energies. To some degree this 
might simply reflect selection rules or structural effects in 
different experiments, but when the same probe sees differ- 
ent states in different experiments (as for 12 C(a,a) 12 C [5, 
6]), there is clearly a problem. It seems safe to infer from 
the data that a broad 2 + -statc exists in the region 14-15 
MeV as e.g. seen in /3-decay [7,1] and in some scatter- 
ing experiments [3], but whether more broad states exist 
near 10 MeV overlapping with the known broad + -state, 
as suggested by [6], is unclear. If so they must be feebly 
fed in the /^-decays as suggested in [11]. The large width 
of the second + -resonance is qualitatively in agreement 
with the results of about 2.7 MeV obtained from the in- 
clastic scattering experiments [5,6]. The existence of nar- 
row states below 10 MeV can be ruled out experimentally. 
The existence of a broad 4 + -state below the known one 
at 14.1 MeV cannot be ruled out experimentally; for com- 
parison the broad + -state centered at 10 MeV is not seen 
in a number of transfer reactions due to its large width — 
it is most clearly seen in /3-decay experiments where a 
4+-state cannot be favorably populated. 

5 Summary and conclusions 

We employ the established method of hyperspherical adi- 
abatic expansion in combination with complex scaling to 
compute the energies and widths of 12 C-resonances below 
15.96 MeV. We use the 3a-cluster model with well known 
angular momentum dependent two-body interactions ad- 
justed to reproduce the low-energy a — a scattering prop- 
erties. Obviously then only three-body properties can be 
computed, but within this overall constraint, the model 
can be applicable even when many-body properties are 
important. The partial waves of each contributing Fad- 
deev component are expanded on a large complete basis 
of Jacobi polynomials. Only very few adiabatic potentials 
are needed for convergence for each angular momentum 
and parity J" . 

The energies are determined by the behavior of the 
potentials at small distances. At intermediate distances 
the potential has a barrier that determines the resonance 
width. As mentioned, the two-body potential is chosen 
to reproduce the low-energy two-body scattering prop- 
erties, and a three-body short-range interaction is also 
added. This accounts for the intrinsic nucleonic degrees 
of freedom by producing a correct boundary condition at 
small distances matching the 3a-structure at large dis- 
tance. However, it does not describe the short-distance 
structure properly and hence electromagnetic transitions 
can not be expected reliably estimated in this model. 
i In our three-body cluster model we add three-body po- 
tentials to place the resonances at the desired positions. 
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These potentials are chosen independently for each J 71 " to 
depend only on the hyperradius. We use one Gaussian 
with a range b = 6 f m corresponding to the distance when 
three a-particles touch each other in an equilateral trian- 
gle. The strength is used to place one of the resonances 
at the desired energy. This invention opens the cluster 
model for applications to many-body resonances similar 
to the simple a-emission model. 

We are not concerned with the bound states where 
especially the ground state can deviate more from a 3a 
cluster state than most of the excited states and here in 
particular from the first + resonance. We focus instead on 
reproducing the correct energy and investigate the system- 
atics of the strongly correlated corresponding width. De- 
signing a three-body potential to fit simultaneously more 
of the states would not provide any more insight. The 
resonance structures are essentially independent of these 
three-body potentials. 

The model leads to a number of low-lying resonances 
for each J n , i.e. two + , three 2 + , two 4 + , and one of each 
of l*, 2~, 3 ± , AT , and 6 + . Dependence on the two and 
three-body interactions are investigated. The conclusions 
are that the highest of the 2 + and 4 + resonances most 
naturally can be placed at the measured values whereas 
the lowest of these resonances acquire a very large width 
when they are pushed towards higher energies. In fact the 
confining barriers of the adiabatic potentials vanish before 
the energy is high enough. 

The structures of the individual resonances are almost 
totally independent of the three-body potential or equiv- 
alently independent of the resonance positions. On the 
other hand, the energies and widths are strongly corre- 
lated and rather insensitive to the structure. With the 
AB(a') Ali-Bodmer potential and the measured energies 
reproduced with a three-body potential of range for b = 
6 fm, the widths are all systematically about a factor of 
two larger than measured indicating a corresponding spec- 
troscopic factor of about 0.5. Two exceptions appear, i.e. 
the second + resonance with a computed width about 
3 times smaller than a recent measurement, and the two 
1+ resonances with 3 — 5 orders of magnitude larger com- 
puted widths and the corresponding very small spectro- 
scopic factors. 

One of each of the negative parity states is found sug- 
gesting that the tentative assignment of one of the 2~- 
resonances should be changed to 4~. The large computed 
width of the second + -resonance appears in the same en- 
ergy region as several other resonances. The two lowest 
2+ and the lowest 4 + resonances are presumably impos- 
sible to reconcile with the experimental information. The 
energies of both the lowest 2 + and 4 + -resonances should 
presumably be higher. This can be achieved by using a 
three-body interaction depending on more than the hyper- 
radius. However, in the present work the focus is on the 
structure and the relation between energies and widths. 

The simple three-body cluster model is, perhaps sur- 
prisingly, doing as well as more microscopic models. The 
difference is really that the phenomenology is inserted on 
the nucleonic level for microscopic models. These micro- 



scopic models are more likely to have problems with spa- 
tially extended systems, and elaborate calculations of ob- 
servablcs. Our cluster model exploits the experimental in- 
formation by inserting the phenomenology on the a clus- 
ter level and a few of the 12 C properties are used as well. 
This model is still technically difficult but closer to the ob- 
servables since we adjust on the three-body level without 
structure changes. 

The discrepancies between experimental and theoret- 
ical results can logically be related to either experiments 
or theory or to a combination of both. Cleaning up on the 
experimental side to achieve consistency between different 
experiments is necessary. This is especially in connection 
with the possibly overlapping resonances around 4 MeV. 
The present model does not treat the intrinsic degrees of 
freedom. They are accounted for on the two-body level 
by phcnomenologically adjusted interactions, and on the 
three-body level by the hyperradial and J 71 ' dependent po- 
tential. The most likely suspect is the simplicity of the 
three-body interaction which is supposed to mock up the 
important effects of the nucleonic degrees of freedom. The 
Pauli exclusion, when three a-particles are close, is per- 
haps not well described. The one-Gaussian structure may 
be too simple, and the potential should perhaps depend 
on angular momentum quantum numbers or on individual 
adiabatic potentials. 

In conclusion, we find in the present model a compa- 
rable number of resonances as in more microscopic calcu- 
lations but the detailed energy spectra vary substantially 
between the different models. We find more low-lying reso- 
nances than measured. The partial wave structures of our 
resonances are very robust independent of their energies. 
We have established resonance structures and suggested 
the most likely energies within the model. 

The real test is then to compare experiment and the- 
ory for observables with strong structure dependence. This 
would demonstrate which of the computed resonances are 
really seen in experiments. The obvious, and perhaps the 
only observable for this is the momentum distributions 
of the a-particles after decay of the resonances. Experi- 
mentally this has to be both accurate and kinematically 
complete. Theoretically large-distance continuum proper- 
ties must be accurately computed and the present tech- 
nique is at the moment the closest to provide answers. 
Both the necessary experiments and calculations are chal- 
lenging but possible and both within reach. 
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Fig. 1. The real parts of the four lowest adiabatic effective potentials, including the three-body potentials, as functions of p for 
the 12 C resonances with J n given in the figures. The two-body interaction, obtained from [17], is a slightly modified version of 
the al-potential of [31]. The parameters (S and b) of the three-body Gaussian potentials, S exp(— p 2 /b 2 ), are given in table 2. 
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Fig. 2. The real (Er) and imaginary (Ei = — -T/2) parts of the energies for ./""-resonances in 12 C obtained after complex 
rotation by the angle 9 also given in the figure. The two and three body interactions are specified in fig. 1. 
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Fig. 3. The real parts of the three lowest adiabatic radial wavefunctions as functions of p for the + , 1 + and 2 + -states of 12 C. 
Their normalization reflects their relative contribution to each resonance. The two and three-body interactions are specified in 
fig. 1. 
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Fig. 6. The partial wave decomposition of the C resonances with .P as indicated on the figure shown as function of p for the 
dominating adiabatic eigenvalue. The two and three-body interactions are specified in fig. 1. 
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Fig. 7. The angular eigenvalues for some of the three-body resonances of 12 C for different interactions. For the + and 2 
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Fig. 8. The real and imaginary parts of the resonance energies for the four cases shown in fig. 7 obtained after rotation by 
the angle 6 for different two-body interactions. The left column is for the two-body interactions "AB(a')" specified in fig. 1. 
The right column uses "AB(a'd')" for + and 2~. The three-body interactions specified by strength and range are given in the 
figures. The lines are the resonance positions followed for b — 6 fm when S is varied 5 MeV up or down from the central value. 
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